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I.  INTRODUCTION 

In  previous  studies,  Dartt,  et  al.  (1974),  Dartt  (1975,  1976),  computer 
programs  have  been  developed  tu  monitor  the  position,  intensity,  and  rota¬ 
tion  of  cyclones  from  digital  satellite  cloud  patterns.  The  current  effort 
provides  a  test  of  these  algorithms  from  a  large  sample  of  tropical  and  extra- 
tropical  storms.  The  technique  is  based  o;v  firJiny  the  distribution  of 
centers  of  curvature  from  elements  of  concentric  cloud  bands  surrounding  the 
storm.  The  location  of  the  maximum  of  this  distribution  is  the  focal  point 
of  the  cloud  bands  and  is  defined  as  the  storm  center.  Other  parameters 
computed  from  this  distribution  yield  objective  measures  of  storm  concen¬ 
tricity  and  rotation. 

A  primary  use  of  the  algorithms  is  to  monitor  cyclonic  vortices  from  the 
high  volume  of  image  data  available  from  geosynchronous  and  orbiting 
satellites.  The  programs  can  be  used  to  follow  the  process  of  storm  develop¬ 
ment  over  the  oceans  v^here  few  conventional  meteorological  observations  are 
available.  In  particular,  any  site  which  has  the  capability  to  receive 
satellite  data  and  has  computer  resources  can  monitor  characteristic  features 
of  storms  In  Its  vicinity. 

Most  emphasis  in  this  investigation  has  been  placed  on  the  analysis  of 
tropical  storms.  Current  operational  techr.iqvies  of  monitoring  tropical  storms 
include  aircraft  reconnaissance  and  satellite  photo  interpretation  by 
experienced  analysts.  The  technique  of  Dvorak  (1975)  of  interpreting  the 
minimum  surface  pressure  of  the  storm  and  the  maximum  wind  from  the  concen¬ 
tric  cloud  pattern.s  is  particularly  useful  for  estimating  extreme  weather 
variations.  Also,  Fett  and  Brand  (1975)  have  analyzed  the  rotation  of 
storms  with  respect  to  changes  in  the  storm's  trajectory..  The  investigation 
conducted  here  Is  to  provide  tools  for  the  support  of  these  efforts.  For 
example,  in  this  report  a  method  is  described  for  obtaining  objective  Dvorak 
T-numbers  from  visual  and  infrared  (IR)  satellite  images  in  order  to  monitor 
the  storm  continuously  over  a  24-hour  period.  Also,  a  reference  axis  is 
defined  from  which  to  measure  storm  rotation  v;ith  time. 

For  extra- tropical  storms,  the  relationships  of  cloud  parameters  are 
analyzed  with  respect  to  position  and  novement  of  the  storm  as  seen  from  a 
satellite.  An  attempt  is  made  to  relate  the  concentricity  o£  the  storm  to 


the  amplitude  of  the  respective  shortwave  trough  at  500  mb.  Nagle  and  Hayden 
(1971)  have  found  good  relationships  between  cloud  parameters  and  shortwave 
troughs  at  this  altitude.  Similarly,  Troup  and  Streton  (1972)  have  classi¬ 
fied  cloud  patterns  of  extra- tropical  storms  in  order  to  initialize  upper- 
air  height  fields  in  data  sparse  regions  of  the  southern  hemisphere. 


II.  CATXULATION  OF  THE  CENTER  OF  CURVATURE  niSTRIBUTION 

A.  TECHNIQUE 

The  procedure  for  computing  the  center  of  curvature  distribution  of 
spiral  cloud  bands  has  been  described  in  detail  in  previous  reports  (Dartt 
c'j  al.,  1974,  and  Dartt,  1975).  It  is  repeated  here  with  emphasis  on  user 
application  to  both  tropical  and  extra-tropical  storms. 

The  first  requirement  is  to  have  an  image  of  a  cyclonic  cloud  pattern 
in  digital  form.  Visual  or  IR  images  can  be  analyzed  using  picture  elements 
of  6-bit  precision  (0-63)  in  gray  scale.  The  image  should  be  scaled  so  that 
the  cloud  bands  have  larger  values  than  the  underlying  sea  surface. 

IR  images  do  not  need  to  be  converted  to  temperature  before  computer  analysis. 

The  basic  image  array  consists  of  128xl28‘*16384  (40000  octal)  elements, 
lliis  is  as  large  an  array  as  is  practical  to  analyze  on  a  general  purpose 
computer  in  batch  processing  mode.  (For  a  128x128  array,  the  program 
storage  is  approximately  140000  octal  locations.)  It  is  important  to  retain 
as  much  image  resolution  as  possible  but  still  contain  the  entire  storm  in 
the  128x128  array.  For  tropical  storms,  an  image  resolution  of  11  Ion  was 
used  in  this  study.  For  extra- tropical  storms,  an  image  resolution  of  22  km 
v:as  u.sed  so  that  the  primary  cloud  vortex  along  with  the  frontal  cloud  band 
could  be  contained  in  the  image  array.  Note;  Changing  the  image  resolution 
results  in  non-linear  changes  in  parameters  derived  from  the  center  of 
curvature  distribution,  and  thus  parameters  derived  2rom  different  resolution 
data  cannot  be  compared. 

For  tropical  storms,  the  tentative  storm  center  of  the  image  should  be 
located  near  the  center  of  the  128x128  array.  Figure  1,  so  that  concentric 
bands  near  the  perimeter  of  the  storm  can  contribute  to  the  center  of 
curvature  distribution.  For  extra- tropical  storms  in  the  northern  hemisphere, 
the  cloud  vortex  should  be  displaced  to  the  left  and  upwards  from  the  center 
of  the  image  array  so  that  the  leading  edge  of  the  frontal  cloud  band  occurs 
within  the  image  array  (see  Figure  17). 

Step  1.  In  the  first  computation  step,  a  digital  filter  is  used  to  enhance 
the  banded  concentric  cloud  structure  of  the  storm.  The  wave  number 

"pon  ;!-  fv’.nctioi!,  S(k  _,k  ),  of  the  filter  is  Gaussian  in  fom  and  radially 
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symnetrical  about  k  •»  k  ■>  0. 
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S(k^,ky)  =  exp 

\  =  Ck^^  +  ky^l^  (1) 

Here  1;  and  k  arc  wave  ritaibers  in  the  x  aiul  v  direction,  k  is  a  comnarabl  ^ 
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radial  x^avc  number,  is  the  radial  wave  number  of  naximum  response,  and  o 
is  a  parameter  controlling  the  width  of  the  response  function.  and  o  are 
chosen  by  the  analyst  to  enhance  the  banded  structure  of  any  desired  scale. 
After  considerable  experimentation  with  both  tropical  (40  images)  and  extra- 
tropical  (11  images)  storms,  the  filter  parameters  Fjn=22  and  o=8  were  chosen 
to  be  optimum  for  producing  distinctive  center  of  curvature  distributions. 

The  peak  filter  response,  F^,  for  11  km  resolution  tropical  images  is  at  a 
wavelength  of  64  km.  For  22  km  resolution  extra- tropical  images,  the  peak 
filter  response  occurs  at  a  wavelength  of  128  km.  The  IR  cloud  pattern  and 
filter  pattern  for  the  western  Pacific  Tropical  Storm  Gilda,  July  04,  1974, 
1226Z,  are  shown  in  Figures  1  and  2,  respectively. 

Step  2.  In  the  second  computation  step,  a  skeleton  of  the  filtered  bands  is 
computed  to  compress  the  filter  pattern  to  a  new  standard  line  segment  form. 
The  orientation  of  the  original  banded  structure  is  preserved  as  seen  in 
Figure  3.  The  edges  and  corners  of  the  filter  pattern  are  discarded  bucause 
of  an  undesirable  halo  sometimes  produced  at  the  image  boundary. 

Step  3.  The  angular  orientation  of  short  line  segments  of  the  skeleton  cloud 
band  pattern  is  now  computed.  This  is  done  by  correlating  a  rectangular  mask 
of  specific  orientation  with  the  skeleton  pattern.  The  mask  is  designed  to 
select  narrow  line  segments  of  fixed  angular  direction  and  inhibit  line 
segments  oriented  in  other  directions.  Repeated  applications  of  this  mask 
for  various  angles  co"\prir.ing  a  360*^  cirolc  g-  ai.rates  a  coMplt^te  angular 
aaal3’sif;  of  thi.-  scg'.ieated  cloud  band  pattern. 

Step  4.  Rays  are  constructed  which  are  orthogonal  to  tlie  directions  of  the 
skeleton  segments  from  the  angular  analysis  above.  The  length  of  the  rays  is 
41  grid  elements  on  either  side  of  the  midpoint  of  the  skeleton  segment. 

For  concentric  cloud  bands,  the  rays  converge  towards  the  center  of  curvature 


and  diverge  on  the  opposite  (convex)  side  of  the  bands.  Ttie  resulting 
pattern  of  converging  rays  from  all  concentric  cloud  band  segments  produces 
a  composite  center  of  curvature  distribution  for  tlie  entire  storm  system. 

This  distribution  is  defined  as  the  number  of  rays  per  cell  within  the  super¬ 
imposed  128x128  cellular  grid.  At  this  stage,  the  computer  has  analyzed  the 
pattern  in  much  the  sn-.e  v:r!y  as  a  human  anilyst,  i.e,,  concentric  cloud  bands 
have  been  selected,  and  the  center  of  curvature  distribution  for  segments  of 
these  bands  has  been  described. 

Step  5.  A  normalization  of  the  center  of  cui'vaturc  distribution  is  now 
performed  so  that  the  final  pattern  primarily  indicates  concentricity  and  is 
not  a  function  of  the  number  of  cloud  band  segments.  If  Steps  2,  3,  and  4 
are  applied  to  a  pattern  of  cloud  bands  with  random  orientation,  the  computed 
center  of  curvature  will  correspond  to  the  center  of  gravity  of  the  pattern. 

To  estimate  the  random  pattern,  a  separate  center  of  curvature  distribution 
is  computed  based  on  the  same  distribution  of  skeleton  segments  as  before. 
However,  a  random  orientation  of  each  segment  is  assumed  rather  than  the 
specific  orientation  found  earlier.  The  resulting  pattern  based  on  the  sum  of 
the  random  components  associated  with  all  skeleton  segments  is  then  subtrac¬ 
ted  from  the  original  center  of  curvature  pattern  to  produce  the  "normalized" 
pattern.  The  normalization  produces  a  positive  pattern  on  the  concave  side 
of  a  cloud  band  and  negative  pattern  on  the  convex  side  of  this  band.  Only 
the  positive  portion  of  the  noimalizcd  center  of  curvature  pattern  is  saved 
for  the  remainder  of  the  analysis.  An  example  of  the  normalized  center  of 
curvature  distribution  for  the  IR  cloud  pattern  of  Gilda  on  July  04,  1974, 
is  shown  in  Figure  4. 

Characteristically,  the  normalization  may  produce  some  noisy  patterns 
in  the  corners  of  the  128x128  grid.  Any  cloud  bands  which  are  concave  with 
respect  to  the  coniuis  '.vil.l  produce  a  reouUial  center  oE  ceivacure 

patrern  tb.ere.  in  Step  ?,  tb.c  cloud  banda  in  the  eornoro  are  rcr.iceed  from 
analysis,  and  the  random  pattern  in  the  corners  is  therefore  always  small. 
Thus,  the  normalized  distribution  of  centers  of  curvature  in  the  grid  corners 
is  larger  than  would  have  occurred  had  the  corner  grid  points  been  retained. 


B.  OBJECTIVE  PARAMETERS 


Various  objective  parameters  of  the  normalized  center  of  curvature 
distribution  are  useful  for  analysis  of  storm  charactcri.stics . 

1.  Locatinf;  the  center  of  the  storm.  Tlie  grid  cell  with  the 
maximum  number  of  intersecting  rays  is  the  focal  point  of  the  concentric 
cloud  hnncls  and  is  defined  as  the  storm  center. 


2 .  MeastU'c  of  stonu  ccncentrfcity  related  to  stomi  intensity . 

The  magnitude  of  t’lc  center  of  curvature  distribution  is  found  by  integrating 
the  distribution  in  a  circular  region  about  the  storm  focal  point  given  above. 
This  magnitude  is  the  average  density  of  rays  throughout  the  circular  area 
and  is  called  the  d  (density) -number  .  Tlie  d-number  is  an  objective  measure 
of  concentricity  of  the  storm  as  indicated  by  curvature  in  the  cloud  bands. 


3.  Reference  axis  for  measuring  storm  rotation.  The  major  axis 
of  the  ellipse  used  to  characterize  the  two  dimensional  center  of  curvature 
distribution  is  a  parameter  indicating  an  axis  of  curvature  for  the  cloud 
bands  which  can  be  monitored  with  tine  from  successive  images.  The  angle  of 
the  major  axis  (y)  with  respect  to  the  coordinate  axes  is  defined  by: 


,  2r  s  X 
..  .  -1  xy  X  V 

2y  =  tan  — ^ - ^2 

s  -  s 
X  y 


(2) 


where  s  and  s  are  the  respective  standard  deviations  of  the  center  of 
X  y 

curvature  distribution  in  the  east-west  and  north-south  directions,  and  r 

xy 

is  the  correlation  of  this  distribution  between  the  x  and  y  directions. 

The  major  axis  is  defined  from  the  same  area  of  the  center  of  curvature 
distribution  surrounding  the  storm  focal  point  used  to  compute  the  d-number. 
The  major  axis  is  indicated  by  a  straight  line  in  the  accompanying  center  of 
curvature  distribution  (Figure  4). 


4.  Stona  el  lipt icily.  To  treasttre  the  rtoim's  elliptical  charac¬ 
teristic,  the  ratio  of  standard  deviations  along  the  major  (o^)  and  minor 
(a.  )  axes  of  the  ellipse  characterizing  the  center  of  curvature  distribution 
is  calculated.  Tliis  parameter  is  useful  for  estimating  the  error  of  position 
of  extra-tropical  storms  as  given  in  IV-C-3. 


III.  TROPICAL  STORM  ANALYSIS 


A.  OBJECTIVE 

Paraiiic:(:ers  of  the  center  of  curvature  distribution  of  cloud  bands  for 
tropical  storms  in  the  western  Pacific  arc  evaluated  for  analyzing  the 
storm's  intensity  (d-menber) ,  rotation  (major  axis),  and  position  (location 

of  riaxin’.um) . 

B.  DATA 

Film  images  of  tropical  storms  were  obtained  from  the  DMSP  satellite 
filra  library  at  the  University  of  Wisconsin,  Madison.  Samples  of  IR  and 
visual  images  of  .selected  typhoons  were  chosen  that  characterized  the  growth 
cycle  of  tropical  storms.  All  storms  selected  had  histories  of  intensity 
and  position  in  the  Annual  Typhoon  Reports  published  by  the  Joint  Typhoon 
Warning  Center  in  Guam.  The  films  were  digitized  with  an  Cptronics  P-1000 
Photoscanner.  \<rliicli  produced  a  basic  500x500  image  array  of  3  km  resolution. 
(The  calibration  of  the  Photoscanner  was  not  changed  during  the  scanning  of 
all  storm  images.)  The  final  128x128  image  array  of  11  km  resolution  used 
for  computation  wa.s  derived  from  this  higher  resolution  array.  Table  I 
shows  the  storm  images  chosen  for  analysis. 

TABLE  1.  Tropical  storm  images  used  for  concentricity  analysis 


Tropical  Storm 

Visual 

Number  of  Images 

IR  Normal  IR 

Inverse 

Nora,  Oct  3-10,  1973 

6 

- 

6 

Patsy,  Oct  7-11,  1973 

4 

- 

7 

Cilda,  Jun  26-Jul  7,  1974 

8 

6 

- 

Mary,  Aug  11-24,  1974 

8 

4 

1 

Polly,  Aug  23-31,  1974 

6 

10 

- 

Uo.;o,  Aug  27-31,  1974 

4 

4 

- 

Shirley,  Sep  4-8,  1974 

3 

4 

2 

The  DMSP  typhoon  images  originally  came  from  variou.s  tactical  sites  in 
the  v.’estern  P.icific  (C'lam,  I’uchu,  Kadcna,  .and  N  akhoii  Phanom)  .  Each  site  has 
its  charactcristi i;  method  of  displaying  iR  data.  In  the  normal  mode,  clouds 


appear  dark  against  a  lighter  background  sea  surface.  In  the  inverse  mode, 
clouds  appear  white  against  a  darker  sea  surface.  All  IR  inverse  mode  images 
except  one  came  from  Fuchu.  For  an  unknown  reason,  the  Fuchu  IR  images  all 
have  a  narrower  range  of  IR  brightness  variability  than  tlic  IR  images  from 
other  sites.  Figure  5  shows  two  IR  image.s  of  Typhoon  Shirley  recorded  within 
a  few  seconds  of  one  another  by  the  same  sensor  (2.4  nautical  mile  [n.rai.]] 
resolution  JK).  The  image  from  Guam  is  o-bit  IK  data  in  the  uomal  mode  but 
has  been  inverted  (brightness  elements  have  been  subtracted  from  63)  so  the 
clouds  appear  white  in  the  usual  meteorological  sense.  The  second  image  is 
6-bit  image  data  fx'Oin  Fuchu,  IR  inverse  mode.  The  graph  below  the  images 
indicates  the  brightness  distribution  of  the  two  images.  Here,  the  bright¬ 
ness  is  also  expressed  as  temperature  using  the  linear  conversion  described 
in  the  DMSP  User  Guide  (1974).  For  September,  the  sea  surface  temperature 
in  this  region  is  approximately  301°K.  According  to  the  D>5SP  User  Guide, 
the  2.4  n.mi.  resolution  IR  sensor  typically  underestimates  the  sea  surface 
temperature  by  5-10°K  in  the  tropics,  depending  on  viewing  angle.  Thus,  the 
image  temperature  distribution  from  Guam  appears  physically  more  reasonable 
than  that  for  Fuchu. 

The  concentricity  analysis  based  on  filtered  cloud  patterns  is  largely 
a  function  of  spatial  variations  in  the  brightness  field  so  differences  in 
calibration  of  the  two  images  do  not  significantly  effect  the  d-number. 
Hov/ever,  in  the  storm  intensity  classification  described  below,  it  is 
necessary  to  define  a  measure  of  area  cloudiness  that  is  dependent  on  an  IR 
brightness  (temperature)  threshold.  Each  type  of  imi-gc  has  a  different 
temperature  range  for  its  gray  scales,  and  different  gray  scale  resolution, 
but  it  should  be  possible  to  compare  absolute  cloud  temperatures.  However, 
in  practice,  this  is  not  the  case  as  explained  later. 

The  IR  image  sample  in  Table  1  contains  a  fev?  images  based  on  the  DMSP 
.33  n.mi.  resolution  sensor  in  addition  to  the  images  based  on  the  2.4  n.ml. 
resolution  sensor.  The  DMSP  User  Guide  indicates  that  the  high  resolution 
sensor  underestimates  temperatures  by  20”k  comp.ared  to  the  low  resolution 
sensor.  Thus,  when  the  .33  n.mi.  data  arc  used  to  determine  area  cloudiness, 
the  20°K  correction  Is  applied. 


All  of  the  above  problems  In  image  and  sensor  calibration  could  be 
avoided  if  a  set  of  precalibrated  digital  data  were  available.  Because  of 
the  high  volume  of  DMSP  data,  a  general  archive  of  all  digital  data  is  not 
feasible.  However,  for  such  phenomena  as  tropical  storms,  such  an  archive 
may  be  practical. 

C.  RliSULTf. 

1.  Tjiph  oon  intensity  c] nsn* f icatioii.  The  Dvorak  technique  of 
classifying  tropical  stoniis  from  satellite  cloud  patterns  is  based  on  an 
analyst's  interpretation  of  storm  concentricity  and  the  area  of  central 
dense  overcast  cloudiness.  The  d-noiuber  defined  from  the  distribution  of 
centers  of  curvature  of  the  storm  is  an  objective  measure  of  concentricity. 
An  investigation  was  undertaken  to  see  if  it  could  be  useful  for  objective 
storm  intensity  classification. 

Initially,  d-numbers  v/ere  computed  from  images  of  three  storms  (Nora, 
Gilda,  and  Mary)  and  compared  to  Dvorak  T-numbers  available  in  the  Annual 
Typhoon  Reports.  These  stonn.a  vjcrc  quite  different  in  their  cloud 
and  Intensity  characteristics.  Nora  was  a  super-storm  with  winds  greater 
than  135  knots  and  a  massive  central  cloud  mass  which  had  few  concentric 
bands  other  than  a  well  defined  eye.  Figure  6a.  Gilda  was  a  tropical  storm 
of  moderate  intensity  and  well  defined  concentricity.  Figure  1.  Mary  was  a 
v/eaker  storm  with  moderate  concentricity.  Figure  6b.  Preliminary  analysis 
of  these  storms  indicated  that  large  aniounts  of  high,  level  clouds  often 
obscurred  the  concentric  banding  of  lov.’cr  level  clouds.  Typically,  in 
the  growth  cycle  of  the  storm,  the  concentricity  as  given  by  the  d-number 
might  be  moderate  (1.0- 2.0)  in  the  early  stages  of  storm  development, 
followed  by  a  decrca.se  (0. 5-1.5)  as  large  amounts  of  high  cloud  were 
generated  by  the  storm.  As  the  storm  reached  its  peak  intensity  and  began 
to  dissip.ite,  tlie  concentricity  increased  ag.ain  and  reached  its  largest 
value  (2. 0-5.0).  'finis,  it  became  apparent  that  an  additional  variable  based 
on  image  cloudiness  \70uld  have  to  be  included  ^?ith  the  d-number  to  account 
for  the  primary  storm  intensity  variability.  Since  both  the  area  of  cloud 
overcast  and  concentricity  are  proportional  to  stom  intensity  in  the  Dvorak 
technique,  it  was  decided  to  forimilate  a  simple  linear  model  of  storm 
intensity  .'ucordii''  to; 


T  =  a  •  (d-munber)  +  b  •  (percentage  image  covered  by  cold  cloud)  (3) 

where  T  Is  the  storm  intensity  as  given  by  the  Dvorak  T-number  and  a  and  b 
are  regression  coefficients  (weights)  to  he  found  from  a  semple  of  d-iroiubers 
and  area  of  cloudiness  calculated  from  satellite  typhoon  images  in  Table  1. 

To  determine  the  statistical  model,  the  percentage  of  the  image  colder 
than  a  thresliolc!  tenipe’.'aturc,  (wliore  T^--220,  230,  ...,  290°K) ,  was 

calculated  for  IR  images  frora  Typh.oons  Cilda  and  Mary.  One  IR  image  from 
Typhoon  Kora  was  also  included  based  on  data  from  Kakhon  Phanon  that  did  not 
indicate  the  calibration  problems  of  the  Fuchu  images.  A  total  of  11  images 
was  used  to  define  the  model.  Separate  models  were  formulated  for  d-numbers 
derived  from  visual  and  IR  images,  although  in  both  models  the  cloudiness 
parameter  is  based  on  the  IR  image.  The  T-numbers  used  for  model  development 
came  from  the  Annual  Typhoon  Report.  For  a  given  satellite  image,  different 
T-numbers  can  occur  because  of  the  interpretation  of  the  cloud  pattern  by 
different  analysts.  The  average  of  individual  T-numbers  was  used  for 
T-number  input  for  model  development  and  also  for  later  model  verification. 

Figure  7  indicates  the  residual  error  of  both  IR  and  visual  models  as 

a  function  of  the  temperature  threshold,  T^.  At  a  threshold  temperature  of 

T  =  270**  (based  on  IR  normal  films) ,  the  error  of  both  models  is  near  minimum 
o 

with  a  stand.'ii'd  deviation  of  approximately  t\;o-thirclc  of  a  T-number.  This 
residual  error  is  larger  than  the  variation  of  T-numbers  by  photo-analysts  as 
described  b>  Sheets  and  Gricman  (1975).  To  reduce  the  error  further  requires 
additional  terms  in  the  model  which  are  not  practical  to  formulate  at  this 
time  because  of  the  limited  data  sample.  For  the  threshold  temperature  of 
270°,  a=.5615  and  b=.0564  for  IR  images;  while  a=.4995  and  b=.0639  for 
visual  images.  (T)ic  units  of  a  arc  T-number/d-niunber  and  for  b  the  units  are 
T-number. ) 

Figure  8  iudici  tj.;  rhe  T-vuuub.-rs  of  the  !;odol  for  Typ;;.,>oa  Clld.’  fvi.-'.u- 
lated  with  2quation  3  and  the  optimized  regression  coefficients  a  and  b 
given  above.  The  individual  contributions  to  the  T-numbers  by  the  concen¬ 
tricity  and  cloud  area  ternts  arc  also  indicated.  The  cloud  area  ajipears 
more  important  than  the  concentricity,  although  the  latter's  importance 
iucrca.sc.s  with  storm  ;:.iturity.  It  can  be  scon  that  hath  the  visu.il  and  TR 


d-numbcrs  jointly  produce  an  internally  consistent  set  of  model  values  whose 
envelope  describes  the  primary  tine  varl<ability  of  observed  T-numbers.  In 
Figure  8,  the  dropsoiide  surface  pressure  from  air  reconnaissance  flights  is 
also  indicated  as  a  direct  measure  of  a  storm  intensity.  Tables  showing  the 
inverse  relationship  hetween  the  minimum  sea  level  pressure  and  satellite 
T-numbers  are  given  in  Dvorak  (1975). 

To  verify  the  model,  the  coefficients  a  and  b  derived  from  the  previous¬ 
ly  described  11  images  are  applied  to  the  IR  and  visual  images  from  Typhoons 
Folly,  Rose,  and  Shirley.  Practically  all  IR  images  from  these  typhoons  are 
in  the  normal  mode  and,  therefore,  consistent  with  the  IR  images  used  to 
define  the  model. 

Polly  was  a  tropical  storm  of  moderate  intensity  (maximum  winds,  95 
knots)  and  good  concentricity  in  its  mature  stage,  Figure  9.  The  model  and 
observed  T-numbers  are  sho^m  in  Figure  10  for  this  storm.  Tlie  envelope  of 
T-numbers  of  the  model  describes  the  primary  storm  variations  quite  well, 
and  there  is  internal  consistency  between  T-numbers  based  on  visual  and  IR 
image  data.  As  in  the  case  of  Gilda,  Figure  8,  the  d-number  contributes 
significantly  to  the  T-number  during  the  storm’s  mature  and  dissipating 
stages.  On  a  systematic  basis,  the  model  appears  to  underestimate  the 
observed  T-niunber  in  the  dissipating  stage  (after  August  29)  by  approximately 
1  to  1.5  T-numbers. 

Tropical  Storns  Rose  and  Shirley  are  distinctly  weaker  stoiras  v;hose 
model  and  observed  T-numbers  are  given  in  Figure  11.  Rose,  Figure  9b,  and 
Polly  were  a  binary  stom  pair  occuring  in  the  same  general  region;  Shirley, 
Figure  5,  occurred  a  few  days  later  in  this  same  area.  The  model  T-numbers 
for  Tropical  Storm  Ro.se  are  v/ithin  approximately  i.5T-nurabers  of  the 
subjective  values.  Ihese  values  for  Shirley  systemacically  underestimate  the 
observed  values  by  about  I  T-number  up  to  the  time  of  maximum  storm  develop¬ 
ment. 

The  verification  statistics  for  Folly,  Rose,  and  Shirley  are  given  in 


Table  2. 


TABLE  2.  Difference  (D),  modeled  minus  observed  T-numbers 


Par.ameter 

IR  Images 

Visible  Tm, 

Average  D 

-.6 

-.4 

RMS  D 

1.0 

.9 

Ln  -  (h)  J  " 

.6 

.8 

Number  of  images 

17 

11 

The  model-derived  T-numbers  appear  too  small  by  approximately  one-half 
T-nuraber  for  both  IR  and  visual  images.  Furthci:,  the  RMS  difference  beUv/een 
modeled  and  observed  T-numbers  is  about  jbl  T-number.  If  the  negative  bias 
between  modeled  and  observed  values  is  removed,  the  resultant  residual  error. 
Table  2-row  3,  is  very  nearly  the  same  as  the  residual  error  of  the  original 
model.  Figure  7. 

There  are  two  sources  of  error  in  the  subjective  analysis  of  T-numbers. 
One  is  the  variation  of  T-numbers  as  a  function  of  analyst,  the  other  is  the 
error  of  approximating  the  minimum  sea  level  pressure  and  maximum  wind  speed 
from  the  T-number.  In  this  investigation,  we  have  not  been  concerned  with 
the  latter  error,  but  rather  with  developing  a  model  that  synthesizes  the 
average  T-munber  of  various  analysts.  Sheets  and  Grlcman  (1975)  indicated 
that  60  to  90  percent  of  T-numbers  produced  by  analysts  from  DMSP  images  are 
within— .5  of  an  official  "best"  T-number.  The  range  reflects  various  analy' 
sis  groups  and  storm  severity  classifications.  Further,  85  to  100  percent 
of  T-numbers  produced  by  analysts  from  this  data  source  are  within  il  of  the 
official  T-nuraber,  and  94  to  100  percent  within  +1.5  of  this  number.  Here, 
10/28  or  36  percent  of  the  comparable  differences  between  modeled  and 
averaged  subjective  T-numbers  lie  within  +.5,  18/28  or  64  percent  within 
+1.0,  and  26/28  or  93  percent  within  +1.5.  Clearly,  the  st:>Listical  model 
described  here  based  on  only  eleven  .storm  iniago.s  (loco  not  have  the  acicurncy 
of  the  human  analyst. 

To  model  T-numbers  from  the  IR  inverse  images  for  Typhoons  Nora  and 
Patsy  (Fuchu  images) ,  a  calibration  of  the  IR  brightness  difference  between 
images  in  the  normal  and  inverse  modes  must  be  devised.  This  was  done  using 
the  tv;o  TR  brightness  distributions  of  the  saviP.  image  of  Typ’u’on  Shirley 


shown  In  Figure  5.  The  threshold  temperature  of  clouds  on  which  the  model 
was  based  was  270°K  for  the  IR  normal  image.  Using  this  threshold  tempera¬ 
ture,  the  percentage  cloud  was  computed  for  the  IR  normal  image  from  Guam. 

For  the  IR  inverse  Fuchu  image,  a  new  threshold  temperature  was  found  which 
gave  the  same  percentage  of  cloud  as  the  previous  image.  This  temperature 
was  238°K  as  indicated  on  the  temperature  distribution  cf  the  IR  inverse 
image.  Figure  5,  The  threshold  temperature  2jS°K  was  then  used  with  the 
rcmaiiiing  images  of  Typhoons  Nora  and  Patsy  to  define  the  percentage  cloudi¬ 
ness  to  be  used  with  the  previous  derived  regression  coefficients. 

Figure  12  shows  the  observed  and  model  T-numbers  for  these  two  storms. 

The  model  underestimates  the  T-number  at  the  time  of  the  storm  maximum  by  as 
much  as  two  T-numbers  for  Patsy,  primarily  because  of  lack  of  concentricity 
of  the  cloud  patterns.  Patsy  and  Nora  were  superstoims  whose  very  cold 
cirrus  shields  were  extensive  and  obscured  the  concentricity  of  lov*er  levels. 
The  DMSP  IR  sensor  does  not  respond  to  temperatures  belov;  210^K,  and  the 
entire  cloud  area  of  these  storms  was  devoid  of  features  in  the  IR  images. 

In  these  storms,  distinctive  diurnal  variations  of  modeled  T-numbers  are 
evident  as  a  result  of  diurnal  variations  in  cloud  area.  The  nighttime 
T-numbers  for  Patsy  (10/07/122,  10/08/122)  and  Nora  (10/05/122,  10/09/122) 
show  distinctive  minima  compared  to  those  computed  at  midday.  Diurnal 
variations  also  trouble  the  subjective  interpretation  of  storm  intensity 
(Dvorak,  1975). 

Because  of  the  tentative  nature  of  calibration  of  these  IR  inverse 
images,  these  cases  were  not  Included  with  the  previous  verification.  The 
average  difference  of  modeled  minus  observed  T-nambers  is  -.7  for  the  IR 
images  and  -.2  for  the  visual  images  indicating  a  systcr:iatic  underestimate 
of  modeled  T-number  as  occurred  in  the  prevtou.*?  verification  based  on  IR 
normal  images.  Only  A/23  or  17  percent  of  the  modeled-observed  differences 
are  within  +.5  T-numbers,  11/23  or  A8  percent  within  +3.0,  and  16/23  or  70 
percent  within  +1.5.  This  degradation  in  model  accuracy  from  the  IR  normal 
images  probably  reflects  the  calibration  problems  between  the  two  film  types. 

Clearly,  at  this  stage  of  the  development,  the  statistical  model  docs 
not  reflect  the  accuracy  of  existing  subjective  techniques.  However,  there  is 
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significant  evidence  of  skill  as  shown  by  the  ability  to  describe  relative 
intensity  changes.  No  doubt  better  results  could  be  obtained  with  a  pre¬ 
calibrated  set  of  digital  data  than  with  the  digital  films  that  were  analyzed 
here.  Further,  there  is  evidence  of  systematic  departures  in  tlie  model  that 
may  be  corrected  with  additional  terms.  Notable  is  the  negative  bias  which 
could  be  corrected  by  a  constant  term.  Also,  a  tendency  to  underestimate 
the  strength  of  stonns  at  maturity  probably  could  bo  corrected  by  including 
a  second  cloud  variable  indicating  the  area  of  high  clouds  or  by  adding  a, 
higher  polynomial  of  concentricity  such  as  the  squa’ e  of  the  d-number, 

2.  Rotation  of  major  axis.  Tne  major  axis  of  the  center  of 
curvature  distribution  is  a  reference  line  that  can  be  mcaitored  for  storm 
rotation,  and  it  is  compared  to  direction  changes  of  the  storm  along  its 
trajectory.  Figure  13  Indicates  the  major  axis  of  curvature  of  visual  and 
IR  images  superimposed  on  the  official  storm  track  of  Typhoon  Mary  (Annual 
Typhoon  Report,  1974).  For  this  storm,  the  major  axes  derived  from  both 
visual  and  IR  images  are  nearly  parallel  to  the  storm  track  for  the  complete 
trajectory  of  the  typhoon.  Figure  14  shows  a  similar  diagram  for  Typhoon 
Gllda.  Here  the  major  axis  is  not  oriented  parallel  to  the  storm  track. 
Still,  changes  in  the  storm's  trajectory  are  reflected  by  the  behavior  of 
Che  major  axis.  For  example,  between  July  1  and  July  4,  the  axis  is  noted 
to  rotate  antlcyclonically  corresponding  to  similar  changes  in  the  storm's 
trajectory.  The  difference  between  the  major  axis  of  visual  and  IR  images 
as  seen  in  Figure  14  is  somewhat  typical  of  most  storms.  The  IR  image 
axis  is  usually  not  the  same  as  the  axis  computed  from  the  visual  image: 
however,  both  axes  generally  show  siniliir  time  behavior  along  the  storm's 
trajectory.  Figvire  15  indicates  the  major  axes  of  the  IR  images  for  Typhoon 
Polly  and  Tropical  Storm  Rose.  Rose  initially  occurred  to  the  west  of  Polly 
and  was  moving  northeastward.  On  August  30,  Rose  was  strongly  influenced  by 
Polly's  circulation  and  began  to  move  rapidly  around  the  southern  flan'.;  of 
this  larger  storm.  Polly  may  .also  have  been  h.!;iuciu:cd  by  lloso  as  Polly's 
trajectory  indlcaces  considerable  cyclonic  curvature  at  this  time.  Ibe  m.ajor 
axis  of  Polly  indicates  consistent  anti  cyclonic  rotation  from  August  24  to 
September  1.  There  are  two  large-scale  anticyclonic  def le.ction.s  in  the  storm 
track  on  August  26  and  August  31.  (Note,  the.  inset  on  Figure  15  indic.ites 


that  the  deflection  on  August  26  may  have  been  cyclonic  V7ith  a  rotation 
greater  than  360°)  Both  direction  changes  occur  when  the  stora's  speed  is 
small  (parameter  A  on  diagram).  The  anticyclonic  rotation  of  the  major  axis 
may  be  a  predictor  of  these  direction  changes.  VJhen  the  storm's  speed  is 
large,  such  as  vhen  it  is  being  .-idvcctod  by  a  preJor.iin int  large-scale  flov7 
pattern,  the  rotation  of  the  major  axis  may  not  be  important. 

Tne  major  axis  analysis  for  Typhoon  Shirley  is  not  sho'w'n  here  but  its 
pattern  V7as  not  unlike  that  for  Typhoon  Gilda  where  rotational  changes  in 
the  axis  appeared  to  occur  concurrently  with  directional  changes  in  the 
storm  trajectory.  Two  storms,  Nora  and  Patsy,  had  very  poor  center  of 
curvature  distributions  because  extensive  amounts  of  cold  clouds  obscured 
the  concentric  banding.  For  these  storms,  the  major  axis  did  not  appear 
related  to  the  storm  track. 

In  conclusion,  the  behavior  of  the  concentricity  axis  of  the  storm 
appears  to  reflect  changes  in  the  storm's  trajectory  at  least  for  storms  cf 
moderate  intensity.  In  one  storm,  Polly,  the  rotation  appeared  to  precede 
the  subsequent  direction  changes.  For  the  remaining  storms,  the  rotation 
changes  in  the  storms'  axis  appeared  to  occur  concurrently  with  similar 
direction  changes  in  the  storm  trajectory  and,  therefore,  did  not  seem  to 
have  any  significant  predictive  value.  As  a  result,  the  major  axis  appears 
most  useful  for  monitoring  current  direction  changes  of  the  storm  and 
perhaps  extrapolating  them  into  the  near  future.  .Mso,  the  major  axis  may 
be  useful  for  helping  to  diagnose  direction  changes  of  storms  where  the 
steering  flow  is  not  well  defined. 

3.  Locating  the  storm  center.  In  an  earlier  study,  Dartt,  et  al 

(1974),  the  location  of  the  maximum  value  in  the  distribution  of  ccf.ters  of 
curvature  was  ompared  to  a  subjective  interpretation  of  the  storm,  center 
for  a  small  number  of  tropical  storm  visual  images.  Here  a  similar  comparis' 
is  made  for  the  large  ntumber  of  visual  and  IR  images  used  in  this  study. 


Table  3. 


TABLE  3.  Cumulative  distribution  of  the  distance  apart 
of  objective  and  subjective  centers  of  tropical  storms 


Image  Number 

IR  A5 

Visual  38 

The  center  can  be  found  more  accurately  from  the  visual  than  the  IR  images. 
This  is  because  the  cloud  tops  of  many  of  the  most  intense  storms  are  so  cold 
as  to  not  exhibit  any  concentric  bands  as  measured  by  the  IR  sensor  on  the 
DMSP  satellite. 

Sheets  and  Grieman  (1975)  discuss  the  analyst's  error  of  storm  position 
from  the  "best  satellite  track".  For  K-ISP  data,  33  percent  of  the  cases 
were  within  20  n.mi.  (37  km)  of  the  optimum  location,  76  percent  were  within 
40  n.mi.  (74  km),  and  91  percent  were  v;ithin  60  n.mi  (111  km).  Thus,  from  a 
comparison  of  the  values  in  Tabic  3  the  subjective  storm  centers  appear 
more  accurate  than  the  objective  centers  found  here.  This  presumes  that  the 
subjective  position  used  for  comparison  with  the  objective  storm  fix  occurs 
along  the  "best  satellite  track".  The  subjective  analyses  are  based  on  the 
very  high  resolution  visual  data  (.6  km)  of  the  DMSP  satellite,  while  the 
objective  technique  uses  data  of  11  km  resolution.  The  accuracy  of  the 
objective  technique  improves  with  higher  rc.se lut ion  data.  In  the  current 
analysis,  the  11  km  resolution  data  v/ere  necessary  for  the  storm  intensity 
model  described  earlier. 

The  difference  between  objective  and  subjective  storm  positions  was 
also  analyzed  a.s  a  function  of  the  d-runr.ber  and  .stPim  cllipticity 

For  storms  v;hosc  d-nu'''.ber  was  less  Chan  1.25,  the  ;ae<ni  position  difference 
was  12  grid  increments  (40  cases).  For  d-nv.rabers  between  1.25  and  2.50, 
the  mean  position  difference  v^as  11  grid  increments  (23  cases),  and  for 
storms  with  d-numbers  greater  than  2,5,  the  position  difference  was  only 
2  grid  increments  (16  cases).  Tlius,  storms  with  good  concentricity  can  be 


Distance  apart  in  11  le.i  grid  units 


0-3  0-6  0-9  0-12  0-15  0-18  0-21  0-24  0-27  0-30  0-33  0-36 


27  44  51  62  73  SO  84  84  89  93  96  1007, 

32  47  63  71  82  84  89  92  95  97  97  lOOli 
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IV.  EXTRA-TROPICAL  STORM  ANALYSIS 

A.  OBJECTIVE 

Paroneters  frotn  the  distribution  of  centers  of  curvature  of  cyclonic 
cloud  patterns  are  analyzed  for  the  purpose  of  calculatil^g  the  position, 
intensity,  and  movement  of  extra- tropical  storms.  The  storms  analyzed  are 
upper  level  vortices  that  typically  have  well-dcf iiicd  spiral  cloud  patterns. 
Such  cloud  vortices  are  well  correlated  vjith  features  of  the  500  mb  geo- 
potential  height  pattern  (Ilaglc  and  Hayden,  1971). 

B.  DATA 

Cloud  concentricity  patterns  are  computed  from  satellite  IR  images  of 
NOAA-4.  Data  for  two  periods,  February  27  -  March  5,  1975,  and  March  26  - 
April  5,  1976,  are  analyzed  for  both  the  Atlantic  and  Pacific  oceans.  The 
procedure  was  to  follow  a  given  storm  through  its  growth  stages  and  compute 
the  respective  parai.ieters  of  the  cloud  patterns.  A  total  of  33  images  for 
five  storms  in  the  Pacific  and  26  images  for  five  storms  in  the  Atlantic 
were  analyzed.  In  the  Pacific,  two  storms  occurred  off  the  east  coast  of 
Japan  and  three  of  the  west  coast  of  California.  In  the  Atlantic,  one  storm 
was  a  weak  subtropical  vortex  and  the  remaining  storms  traversed  the  ocean 
at  higher  latitudes. 

Figures  16  and  17  show  one  of  the  extra- tropical  storms  in  the  western 
Pacific  at  successive  12-hour  intervals.  The  distribution  of  the  centers  of 
curvature  evolve  in  much  the  same  way  as  for  a  tropical  storm,  and  for  this 
storm  the  intensity  given  by  the  d-number  attains  the  large  values  that  are 
characteristic  of  the  tropics.  Typically  the  center  of  curvature  distribu¬ 
tion  for  extra-tropical  storms  is  more  elongated  than  circular  as  a  result  of 
the  comma-shaped  cloud  pattern  associated  with  the  frontal  band. 

C.  RHSULTS 

1.  Comparison  of  d-number  and  500  mb  trough  magnitude.  The 
hypothesis  was  made  that  concentricity  in  the  cloud  patterns  as  measured 
by  the  d-number  is  correlated  to  geostrophic  vorticity  computed  Irom  the 

geopotcntlal  height  field  (z); 
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f  =  2Sl  Sint-*  (4) 

where  g  in  the  gravitational  constant  and  f  is  the  Coriolis  parameter 
derived  fro:a  the  roLaCion ’r ate  oC  the  earth.  end  latitude  (G*) .  'Xlie 
Laplacian  of  a  scalar  field  is  a  measure  of  the  departure  of  a  local  value 
from  the  surrounding  regional  average.  Thus,  the  gcostropViic  vorticity  is 
well  correlated  with  the  shortv.’ave  portion  of  the  original  height  field. 

As  a  result,  the  product  of  the  d-numher  and  sinC*  should  be  correlated  v/ith 
the  amplitude  of  the  shortv/ave  trough.  The  amplitude  of  the  shortwave  trough 
was  computed  by  application  of  the  Moll  (1963)  scale  analysis  program  to 
KEPRF  grid  data  for  12-hour  periods  corresponding  to  the  intervals  of  - 
satellite  images.  This  program  partitions  the  variability  of  the  height 
field  into  long  and  short  wave  components.  The  spectral  characteristics  of 
the  shortwave  component  are  controlled  by  a  parameter  or.  For  this  study 
a=5,  and  the  amplitude  of  waves  shorter  than  five  grid  lengths  (one  grid 
length  is  361  km)  are  unchanged  by  the  analysis.  However,  waves  larger  than 
24  grid  lengths  arc  completely  removed,  and  the  amplitude  of  waves  of  ll  grid 
lengths  is  reduced  by  50%. 

As  a  first  test,  five  possible  cloud  predictor  variables  were  correlated 
with  the  amplitude  of  the  500  mb  shortwave  trough.  These  included  latitude, 
©,  cloud  diameter,  D,  cloud  amplitude.  A,  concentricity,  d* sinO  ,  and 
percentage  area  cloudiness,  C,  where  C  is  defined  as  the  fraction  of  the  IR 
temperature  pattern  colder  than  a  given  threshold  Tq.  Nagle  and  Hayden 
(1971)  found  the  first  three  variables  to  have  significant  correlation  with 
the  500  mb  shortv/ave  pattern  over  a  much  larger  data  sample  than  described 
here.  The  cloud  diameter,  D,  corresponds  to  an  east‘\7est  dimension  of  the 
cyclonic  cloud  pattern,  while  ti;o  amplitude,  is  a  dii.ienslcn  from  the  cloud 
vortex  south  to  the  frontal  band.  Both  D  and  A  are  estimated  from  satellite 
cloud  images  in  Environmental  Satellite  Imagery,  using  procedures  described 
by  the  authors.  These  variables  were  included  as  standards  to  which  corre¬ 
lations  of  other  variables  could  be  compared  for  this  smaller  data  sample. 

All  correlations  were  stratified  by  wliethcr  the  cyclones  occurred  in  the 


Atlantic  or  Pacific  ocean.  Table  4  Indicates  correlations  for  the  years 
1975  and  1976  separately. 


TABLE  4.  Correlation  of  cloud  variables 


and  the  cimplitucie 

(2)  of 

tl.e  500 

l  ib  .shor 

tv.v’.vt:  trouj 

Data 

Z«Lat 

ZMJ 

Z-vA  ’L 

•d  sin'"' 

Z»C  (T  =240^10 

Atlantic  (1975-14  images) 

.89 

.33 

-47 

.63 

•  O 

.63 

Atlantic  (1976-11  images) 

-.50 

.50 

-.41 

-.18 

.72 

Pacific  (1975-21  Images) 

-.01 

.08 

.38 

-.32 

.11 

Pacific  (1976-12  images) 

.10 

.22 

.32 

.02 

-.31 

The  regional  and  yearly  differences  are  substantial.  Initially,  it  was 
intended  to  use  relationships  found  from  correlations  in  1975  to  develop  a 
model  for  specifying  shortwave  trough  amplitude.  The  model  would  then  be 
verified  using  cloud  patterns  from  1976.  llov/ever,  the  regional  and  yearly 
differences  are  so  marked  that  the  formulation  of  a  general  linear  model 
was  not  possible. 

Differences  in  correlation  between  variables  can  occu.r  because  of 
phase  differences  of  these  variable.-,  v/lth  time.  Since  most  of  our  images 
are  time  sequences  characterizing  the  growth  cycles  of  individual  storms, 
it  is  possible  to  an'ilyzc  dilfereuccs  in  the  phase  of  the  two  variables. 
Figure  18  is  a  time  composite  of  d'siuo  and  C  w’ith  respect  to  the  time  of 
maximum  trough  amplitude  Iti  the  Atlantic,  the  concentricity  as  given 

by  the  d-number  is  in  phase  v/ith  the  trough  minima;  but  in  the  Pacific,  the 
maximian  d-number  occurs  approximately  24  Ijo’.irs  prior  to  the  trough  mluimua. 
In  the  Pacific,  for  each  of  tlie  five  .storp'.  sequencfs  apalyzed,  the  maxinum 
d-number  occurs  prior  to  maximum  trough  development.  In  the  Atlantic,  the 
maximum  concentricity  occurs  first  only  two  out  of  five  times.  If  concen¬ 
tricity  is  a  precursor  of  trough  development  in  the  Pacific,  then  this 
variable  might  be  useful  for  prognostic  purposes.  However,  it  docs  not 
seem  to  be  useful  by  itself  as  a  single  vnri.-ible  for  classifying  a  storm  on 


An  alternative  explanation  of  differences  in  the  correlations  may  be 
caused  by  the  difference  in  observational  coverage  bet.wecn  the  Atlantic  and 
Pacific  oceans.  The  growth  cycle  of  Pacific  stoms  occurs  in  regions  of  poor 
data  coverage  and  maturity  occurs  near  land  where  upper  air  data  are  present, 
la  regions  of  missing  data,  nrost  tecliniqncs  of  data  interpc>lat J on  are  C'.iii- 
serv'ative  v;hen  specifying  the  initi.al  fields.  Thus,  in  remote  .sreas,  the 
araplitude  of  troughs  is  likely  to  be  underestimated.  If  this  is  the  case, 
then  the  NEI'RF  height  field  used  here  for  correlation  is  likely  to  reflect 
the  same  systematic  behavior.  Under*  these  circuiustanccs ,  the  use  of  objec¬ 
tive  concentricity  to  help  specify  in-phase  trough  development  would  be  very 
useful  indeed. 

2.  Relation  of  major  axis  to  storr.i  movement.  Figures  19  and  20 
show  the  major  axis  of  concentricity  for  four  storm  sequences  in  the  Pacific 
and  two  storm  sequences  in  the  Atlantic.  Included  on  these  di.'igrnnis  is  the 
trough  axis  Indicated  by  a  subjective  analysis  of  curvature  in  the  cloud 
bands  to  the  south  of  the  vortex  center.  Much  of  the  time,  there  appears 
to  be  a  correspondence  between  the  axis  of  the  trough  and  the  major  axis. 

Other  times,  the  axes  are  systcmc4tically  displaced  from  one  another  by  a 
fixed  amount.  The  axis  of  concentricity  primarily  reflects  the  tight  curva¬ 
ture  of  the  cloud  bands  to  the  north  of  the  cloud  vortex,  and,  therefore, 
need  not  correspond  to  the  trough  axis.  Table  5  indicates  the  preference 
of  storm  movement  as  a  function  of  orieiitation  of  the  major  axis.  Wien  the 
orientation  of  the  axis  is  between  150-180°  and  0-30°,  there  ii;  a  high 
probability  that  the  storm  will  move  cast  or  southeast.  If  the  major  axis 
is  v;ithin  60°  of  east  (30-150°),  the  storm  will  very  likely  move  to  the 
northeast. 

TAI'I.E  5.  Distribution  of  t\.'clvc-hour  storm  t:  -'.j'.c.  tosi  c-.. 

as  a  function  of  preceding  orientation  of  the  major  axis 

Major  Axis  Orientation 

Storm 

Movement  180  to  150°  150  to  120°  120  to  90°  80  to  60°  60  ro  30°  30  to  0 


Northeast  2  5  2152 


.Southeast 


3.  Locating  the  storm  center.  Extra- tropical  storms  are  less 
circular  than  tropical  storms  and  have  a  less  distinctive  center  of  curva¬ 
ture  di£;tribution,  Iho  ratio  of  the-  standard  deviations  of  variability 
along  the  axes  characterizing  the  center  of  the  curvature  distribution  is 
a  measure  of  the  ellipticity,  Il-h  4,  and  the  d-number  is  an  objective  measure 
of  concentricity,  both  parci'v.eters  arc  useful  for  classifying  the  error  in 
C'bject'ivi  ly  locating  the  center  of  extra- tropical  .'^torris.  The  average 
difference  be  Lvrecn  th.c  subjc.ctive  and  objc-ctivc  storia  centers  is  given  in 
Table  6  as  a  function  of  these  tv;o  variables 

TABLE  6.  Objective  and  subjective  storm  center  differences  (x22  km) 
(Numbers  in  parenthesis  are  the  number  of  images  examined) 


Ellipticity  (a  A’.  ) 

cl  D 


Q-nuniDor 


0-1.25  r 1.26-2.50  i  2.50 


1.0-1. 5  9(18)  8(8)  1(4) 

1. 5-2.0  15(17)  8(7) 


A  relatively  poor  center  (15x22=330  km)  is  found  when  the  d-number  is 
small  and  the  ellipticity  is  high.  Conversely,  a  very  good  center  is  found 
when  the  d-n\.unber  is  large  and  the  storm  is  nearly  circular  (1x22=22  km). 
Increasing  the  resolution  of  the  basic  cloud  pattern  image  data  would 
decrease  these  errors. 
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SUMMARY  AND  CONCLUSIONS 


Techniques  that  have  previously  been  developed^^for  monitoring  the 
intensity,  location,  and  rotation  of  cyclonic  storms  from  digital  satellite 
cloud  patterns  are  tested  in  this  study .,for  a  sample  of  83  tropical  and  59 
extra-tropical  cyclone  images.  The  caiputer  programs  are  designed  to  produce 
a  nonaalized  dist  rih'j.tioii"  of  cc.itors  oi;  curvature  for  coticenUric  clouds  bands 
surrounding  a  storm  Tl'.c  location  of  the  maximum  of  this  distribution  is 
used  as  a  measure  of  storm  position.  Tiic  horizontally  integrated  magnitude 
of  the  distribution^  tbc  d-niimher^  is  reed  as  a  measure  of  storm  intensity, 
and  the  major  axis  of  th.c  distribution  is  used  as  a  reference  line  to 
measure  storm  rotation. 

A.  TROPICAL  STORliS 

For  western  Pacific  typhoons,  the  cyclone  image  d-number  was  used  with 
an  area  cloudiness  parameter  to  formulate  an  objective  intensity  parameter 
comparable  to  the  Dvorak  T- number.  The  T-uu:abor  is  currently  used  opera¬ 
tionally  to  describe  storm  intensity  from  satellite  cloud  patterns.  A 
portion  of  the  sample  images  was  used  to  construct  a  statistical  model  of 
T-nuinber  variability  using  parameters  computed  from  vi.sual  and  IR  images. 

The  model  appears  to  describe  the  primary  tine  variations  of  intensity 
charactcrij'.ing  the  growth  cycle  of  tropical  storms.  In  terms  of  absolute 
accuracy,  the  difference  betv;eon  the  model  and  subjective  T-numbers  is 
approximately  +1.0.  Corresponding  accuracy  of  T-numbers  by  analysts  is 
approximately  +.5.  To  improve  the.  statistical  model,  recommendations  are 
made  to  add  terms  to  account  for  systematic  differences  between  the  model 
and  observed  T-numbers. 

The  rotation  of  the  major  axis  of  concentricity  of  the  storm  appears  to 
be  correlated  with  direction  changes  of  the  storm's  trajectory.  Generally, 
anCicyclonic  (cyclonic)  cl-.angco  in  rotation  of  tiie  major  axis  occur  concur¬ 
rently  v/ith  anticycJoaic  (cyclonic)  changes  in  .storm  direction.  For  one 
typhoon,  Polly,  it  appeared  that  the  rotation  n.ay  have  preceded  the  subse¬ 
quent  direction  changes  of  the  storm's  trajectory.  The  major  axis  may  be 
useful  for  diagnosing  the  movement  of  erratic  storms  where  a  large-scale 
steering  flow  is  imt  well  clcfin.ad. 


The  objective  tropical  storm  center  compares  favorably  to  the  storm 
position  specified  by  the  analyst  for  storms  of  strong  concentricity.  For 
storms  vjith  few  concentric  bands,  the  analyst  can  position  the  storm  more 
accurately.  An  improvement  in  the  resolution  of  the  satellite  data 
analyzed  hero  (11  km)  v/ill  improve  the  accuracy  of  the  obiective  tr-ehniene. 

B.  E>:TRi\-TllO?j;CAL  STOil'iS 

The  relationship  betv’oen  concentric  storm  propei-tics  of  oxtra-trop Leal 
storms  and  the  500  irb  height  pattcr-n  was  much  le.ss  distinct.  Regional  and 
yearly  differences  in  correlations  betv/een  cloud  parameters  and  the  ampli¬ 
tude  of  the  500  mb  shortwave  trough  were  substantial  and  prevented  any  model 
development. 

In  a  diagnostic  investigation  of  the  regional  differences  of  correlations 
it  was  found  that  for  storms  in  the  Pacific  the  maximum  concentricity  (d- 
number)  precedes  the  development  of  the  trough.  In  the  Atlantic,  the  tendency 
was  for  the  maximum  d-number  to  occur  simultaneously  with  trough  development. 

The  major  axis  of  concentricity  of  extra-tropical  cloud  patterns  appears 
to  generally  reflect  the  curvature  of  the  cloud  bands  to  the  north  of  the 
cloud  vortex.  Often  this  axis  is  coincident  with  the  axis  of  the  trough  as 
seen  by  a  s^ihjcctive  analysis  of  curvature  in  the  cloud  bauds  to  the  south  of 
the  cloud  vortex.  A  tendency  is  noted  for  storms  to  move  .southeast  or  east 
if  the  major  axis  is  within  +30°  of  south.  If  the  major  axis  is  oriented 
^0°  of  east,  the  storm  v/ill  very  likely  move  norVheastwai d  in  the  next  12- 
hour  period. 

Finally,  an  evaluation  of  objective  storm  position  was  made  for  extra- 
tropical  storms.  The  accuracy  of  cloud  vortex  position  is  dependent  on  the 
clrcu.larity  of  the  storr.’  and  the  d-menber.  Extra- tropical  storms  r.ca!o!-  be 
objectively  posit  i.nnod  ns  .accurately  an  fopical  storiu.:. 

C.  RCCroLMENDATlO.NS 

The  most  promising  application  of^ the  objective  techniques  is  lor  ' 

the  analysis  ei^the  intensity  and  rotation  of  tropical  storms.  The  results 
described  in  this  report  arc  based  on  UMSP  film  images  of  tropical  storms  ■  y 


which  were^thcn  digitized.  There  are  problems  in  calibrating  the  various 
types  of  -Ht,  images  in  terms  of  temperature.  The  regression  coefficients  used 
in  the  storm  intensity  model  and  the  threshold  v.tIuc  of  temperature  used  in 
defining  percentage  cloudiness  may  not  be  applicable  to  ap.other  data  set. 

A  technique  to  calibrate  IR  brightness  between  satellites,  perhaps  based 
on  sea  surface  tenperatur?,  should  be  developed.^  Any  operational  applica¬ 
tion  of  the  above  technique  will  bo  performed  on  lprecali:>r<itc.d  digital,  image 
data.  Tnorefore,  add:' tiotml  model  doveloiraent  aiicl  subsequant  vf;ri.fication 
should  be  performed  on  such  precalibxa*  ed  digital  data  be  .'oris  eventual 
application.  The  techniques  developed  here  vjith  photographic  data  should  be 
verified  with  direct  digital  data. 
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Ficure  2.  Typhoon  Gllda  IR  filter  pattern  04  July,  1226Z 
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Figure  3.  Typhoon  Gilda  lU  skeleton  pattern  O.'*  July,  1226Z 
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Figure  4.  Typhoon  Gilda  IR  center  of  curvature  pattern  0^  July,  1226E 


Nuisber  brightness  elements 


Guam  inverted  IR  normal 


Fuclm  Hi  Inverse 


4  a 


Guam  ,’i 

V.../ 


'  >1 


\ 

« 

V 

■  "I 

60 


296  K 


Figure  5.  DMSP  IR  images  for  Typhoon  Shirley,  Sep  6,  1974  (12102)  from 

Guam  and  Fuchu  sites  and  the  respective  IR  brightness  distribu¬ 
tion  of  each  image. 
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Figure  7.  Residual  error  of  statistical.  r.iod''l  as  a  function  of  temperature 
threshold,  ’i’  ,  for  defining  percentage  cloudiness.  Uot.s  indicate 
the  error  c«irve  for  visual  images,  circles  the  error  £or  IR 
imuge.s . 


number 


Figure  8.  Intensity  history  of  Typhoon  Gil  da  as  given  by  (a)  observed  .  - 
numbers  (dots)  and  dropsonde  surface  pressure  (circles)  and  by 
(b)  T-nuitibcrs  calcv;latecl  by  the  r-tu<lel  from  visual  (solid  line) 
and  IR  (dashed  lino.)  images.  Tlic  shaded  area  is  the  contribu¬ 
tion  to  the:  T-nVih.'.-r  in  wud;.!  from  the  percont  cloudiness 
v.ir  j  able. 
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Figure  10.  Intensity  history  o£  Typhoon  Polly  as  given  by  (a)  observed  T- 
numbers  (dots)  and  dropsonde  surfr-co  pressure  (circles)  and  by 
(b)  T  -numbers  calculated  by  the  model  from  visual  (solid  line) 
and  IR  (dashed  line)  images.  The  shaded  area  is  the  contribution 
to  the  T-number  of  the  model  from  the  percent  cloudiness  variable 
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Figure  11.  Intensity  history  of  Typhoon  Shirlc^y  a'vl  Trop ! ral  SLoip  t; 
■  given  by  (a)  observed  T-numbevs  (dot.s)  and  Oropaoin'o  aurC 

pressure  (circles)  and  by  (b)  T-uur.!b«  i.s  cnlcvi .'.Led  by  lln 
from  visual  (solid  line)  and  IK  (daslud  lin-,)  iwages.  Sh 
.’I'ca  is  ti'.e  eontrl’iut  i  I -I  Vo  Lb.  -  T- r.v;  ■  '  o,*'  ;'.■■■  ;  i.^ 

I'Crccnt  cloudiness  variable. 
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Figure  12.  Intensity  histories  of  Typhoons  Patsy  and  Nora  as  given  by  (a) 

observed  I’-nuiTibcrs  (dots)  and  dropsonde  surface  pressure  (circles) 
and  by  (b)  T-nunbers  calculated  by  the  model  from  visual  (solid 
line)  and  JP.  (dashed  lino)  images.  Shaded  area  is  the  contribu¬ 
tion  to  the;  T-ne -.her  of  the  model  by  the  percent  cloudiness 
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Figure  16.  Infrared  Images  and  center  of  curvature  distributions  for  extra- 
tropical  storm  In  the  western  Pacific  on  (a)  March  27,  1976  (12Z) 
and  (b)  March  28,  1976  (OOZ). 
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Figure  17.  Infrared  Images  and  center  of  curvature  distributions  for  extra- 
tropical  storm  in  western  Pacific  on  (a)  March  28,  1976  (IIZ)  and 
(b)  March  28,  1976  (237.). 
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Figure  18,  Time  relation  of  <  oncentvic  Lty  (d  r.in  0)  and  cloudiness  to  the 

minimum  of  the  500  mb  shertvave  trough.  Dots  indicate  variables 
for  storms  in  tho  Pacific;  circles  indicate  variables  for  .storm.s 
in  the  Atlantic.  Tlie  time  of  trough  mini!.ium  is  t„ 


the  major  axis 
trough  (dashed 
ring  1975  and  19 


or  of 
level 
ic  du 


